Chromoelectric oscillations in a dynamically evolving anisotropic background 
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I. INTRODUCTION 

The purpose of ongoing and upcoming heavy ion col- 
lision experiments at the Relativistic Heavy Ion Collider 
(RHIC) and the Large Hadron Collider (LHC) is to study 
the behavior of nuclear matter at high energy density, 
e > 1 GcY/ira. At such high energy densities one ex- 
pects to create a deconfincd quark gluon plasma (QGP). 
With such experiments one hopes to not only cross the 
threshold necessary to create a QGP, but to also study 
its properties such as transport coefficients, color opacity, 
etc. One complicating factor is that the QGP generated 
in such collisions lasts for only a few fm/c and during 
this time the bulk properties of the system, e.g. energy 
density and pressure, can change rapidly. Therefore, dy- 
namical models that can describe the evolution of the 
system on the fm/c timescale are necessary in order to 
make reliable phenomenological predictions. 

To first approximation, it seems that the dynamics of 
the soft background is well-described by relativistic vis- 
cous hydrodynamics P^13j . However, viscous corrections 
to the ideal energy momentum tensor cause it to become 
anisotropic in the local rest frame of the system. For 
small deviations from isotropy, 2nd-order viscous hydro- 
dynamics describes the evolution quite well; however, for 
large deviations from isotropy this is no longer the case. 
Large deviations from isotropy occur at very early times 
after the initial nuclear impact and near the transverse 
or longitudinal edges of the plasma where the matter 
is nearly free streaming. The presence of momentum- 
space anisotropies seems unavoidable in dynamical mod- 



els. In fact, even in the limit of infinite strong coupling, 
momentum-space anisotropies persist during the entire 
lifetime of the plasma [l4lI16j . Large momentum-space 
anisotropies pose a problem for 2nd-order viscous hy- 
drodynamics since it relies on a linearization around an 
isotropic background. If the linear corrections grow too 
large this can generate unphysical results such as nega- 
tive particle pressures, negative one-particle distribution 
functions, etc. [T7] . 

In order to ameliorate these problems it is possible 
to reorganize the derivation of the necessary dynami- 
cal equations by linearizing around an anisotropic in- 
stead of isotropic background. Doing so results in a 
dynamical framework called anisotropic hydrodynamics 
[I8l424j . In the limit of small deviations from isotropy, 
anisotropic hydrodynamics reduces to 2nd-order viscous 
hydrodynamics, but can also faithfully describe large de- 
viations from isotropy such as those created during the 
initial longitudinal free streaming phase of the plasma 
lifetime. This framework has now been used to model the 
full (3-|-l)-dimensional dynamics of the QGP |24j. Com- 
parison of the anisotropic hydrodynamics predictions for 
observables such as the bulk flow as a function of trans- 
verse momentum and rapidity with experimental data 
indicate that it is possible that large momentum-space 
anisotropies can persist for up to 2 - 3 fm/c after the 
initial nuclear impact. Given this, it is imperative to 
revisit the study of basic properties of the QGP in a 
time-evolving anisotropic background. 

In this paper we study the oscillations of a uniform lon- 
gitudinal chromoelectric field in a dynamically-evolving 
momentum-space anisotropic background in the weak 
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field limit. For simplicity, in this work we restrict 
ourselves to a (O-l-l)-dimensional boost-invariant back- 
ground. The necessary anisotropic hydrodynamics equa- 
tions are obtained from the first two moments of the 
Boltzmann-Vlasov equation using a spheroidal form for 
the one-particle distribution function in the local rest 
frame |25| . The dynamical equations in this case were 
first obtained in Refs. [HI In both Ref. [Hj and 

Ref. [inj a timescale for the approach to isotropic ther- 
mal equilibrium, Teq, was introduced. In Ref. [18] r^q 
was assumed to be a constant, while in Ref. [12] this 
time scale was proportional to the average local inverse 
transverse momentum of the plasma constituents and was 
determined self-consistently in terms of the local plasma 
environment. 

In the case that Tcq is constant, the late time behav- 
ior of the system is that of ideal hydrodynamics. In the 
case that r^q is proportional to the local inverse trans- 
verse momentum scale, the proportionality constant can 
be fixed by requiring that the late time dynamics of 
the system is that of 2nd-order viscous hydrodynamics. 
We will consider both cases in order to assess the im- 
pact of this choice. In either case the anisotropic hy- 
drodynamics equations provide the proper-time depen- 
dence of the local transverse momentum scale, A(r), 
and momentum-space anisotropy ^(r) = — 1 

where (p^) and (py) are the average transverse and lon- 
gitudinal (beamline-direction) momenta squared in the 
local rest frame of the plasma constituents, respectively. 

Given this time evolving background, we linearize the 
Boltzmann-Vlasov equation in order to study the evolu- 
tion of a uniform longitudinal chromoelectric field fiuc- 
tuation, £z. We consider the weak-field limit in which 
case we can use the abelian dominance approximation 
for the color fields [26 29' (see also Ref. [3D]). In a 
static constant-temperature plasma, uniform longitudi- 
nal field fluctuations oscillate in time with a frequency 
given by the plasma frequency Wpi — niD/ VS where 
mjj = {Nc/3 + Nf /6)g'^T^ is the leading-order gluonic 
Debye mass. In a time evolving system, the plasma fre- 
quency is time dependent and one must self-consistently 
solve the linearized Boltzmann-Vlasov equation together 
with the Maxwell equations. In general, the result can 
be cast in the form of an integro-differential equation for 
the evolution of Sz- For the case of ideal hydrodynamical 
evolution, Bialas and Czyz |31| derived such an equation 
and solved it numerically. 

Here we extend the treatment of Bialas and Czyz 
to (i) include a dynamically evolving anisotropic back- 
ground and (ii) include the effect of collisional damp- 
ing. We will present numerical solutions to the resulting 
integro-differential equations for small and large magni- 
tude momentum-space anisotropics in order to assess the 
impact of momentum-space anisotropy on plasma oscilla- 
tions. The equations we obtain are applicable to an arbi- 
trary time-dependent anisotropic background. Although 
we consider the evolution of a stable longitudinal chro- 
moelectric field, the techniques used herein could have 



application to the study of the evolution of non- Abelian 
plasma instabilities in a dynamically evolving anisotropic 
background [35H31]. As a cross check of our results we 
present a comparison with the results of Ref. which 
presented an analysis of all stable and unstable collective 
modes of the QGP in the limit of a longitudinally free 
streaming background. We show that in this limit we 
obtain the same evolution and asymptotic behavior as 
in Ref. [SI], giving us confidence in our theoretical and 
numerical methods. 

The structure of the paper is as follows. In Section 
|TT|we specify the conventions we will use throughout the 
paper. In Section III we review the semi-classical trans- 



port equations for the quark gluon plasma in the abelian 
dominance approximation. In Section |IV| we linearize 
the Boltzmann-Vlasov equation to zeroth and first or- 
der in fiuctuations. In Section [Vj we couple the fiuctua- 
tions via currents to the Maxwell equations and solve the 
coupled Boltzmann-Vlasov-Maxwell system of equations 
to obtain an integro-differential equation which governs 
the time evolution of uniform longitudinal chromoelec- 
tric fluctuations. In Section \VJ\ we present the results of 
numerical solution to the evolution equations for differ- 
ent types of anisotropic backgrounds and compare with 
numerical and analytic results available in the limit of 
longitudinal free streaming. In Section |VII| we present 
our conclusions and an outlook for the future. 



II. CONVENTIONS 

Below we use the following deflnitions for momentum 
rapidity (y) and spacetime rapidity (77), 



1, ^ + 

w = - In 

2 E-pn 



I, t + z 

77 = - In , 

' 2 t-z' 



(1) 



which come from the standard parameterization of the 
four-momentum and spacetime coordinates of a particle, 

pf^ ^ (E,pj_,p\\) = (m_L cosh y,p±,mj_ sinh y) , 

= {t,x±, z) — (r cosh r],x±,T sinh 77) . (2) 

In Eq. ([2]) the quantity mj^ is the transverse mass 

m± = ^Jm'^ +pI+pIi (3) 

and r is the proper time 

T = -z2. (4) 

Throughout the paper we use natural units where c — \ 
and fi — 1. 



III. SEMI-CLASSICAL KINETIC EQUATIONS 
FOR QUARK-GLUON PLASMA 

In the abelian dominance approximation, the transport 
equations for quarks, antiquarks, and gluons have the 



3 



form 

{p^d^ ± ge, • F'^'-p^dP) Qt {x,p) = Cf , (5) 

{P'd,. + 5^?,, • F^^'Pudl) G., {x,p) = C., , (6) 

where Qt{x,p), Q^{x,p), and Gij{x,p) are the phase- 
space densities of quarks, antiquarks, and charged gluons, 
respectively. Here g is the strong couphng constant, and 
i,j = (1,2,3) are color indices. The terms on the left- 
hand-side describe the free motion of the particles and 
the interaction of the particles with the mean field F^^,. 
The latter describes neutral gluons [35] . 

In this work, the only non-zero components of the ten- 
sor F^'' — (-P'(3) , ) are those corresponding to the 
longitudinal chromoelectric field £z — {F^^y F^^^^). The 
quarks couple to the chromoelectric field £z through the 
charges 




The gluons couple to €z through the charges rj^j defined 
by the relation 

Vij=^t-^r (8) 
Below, we make use of the following relations 

Ea b ^ rah 
^^^^ = 2^ ' 



(9) 



where a, 6 G {3, 8}. 

The terms on the right-hand-side of Eqs. ^ and ^ 
are collisions terms, which we treat in the relaxation time 
approximation 



Cf=-p'^U^- 



'-{x,p) - QtJx^p) 



r< _ _„prr G^jix,p) - Gcq(a;,p) 



(10) 
(11) 



' cq 



Here is the four-velocity of the local rest frame 

U^^-f{l,v^,Vy,Vz)n={^-v^r^- (12) 

In this paper we consider boost-invariant longitudinal ex- 
pansion and hence we set Vx = Vy = Q and Vz = z/t. 

IV. LINEARIZATION OF KINETIC 
EQUATIONS 

In the following, we seek the solutions of kinetic equa- 
tions of the form 

Qf{x,p) = Q^{x,p) + 5Qf{x,p), (13) 
Gy(a;,p) = Go{x,p) + SG^j{x,p), (14) 



where the corrections to the background distributions are 
proportional to the coupling. We emphasize that the 
background distributions Qq{x,p) and Go{x,p) are dif- 
ferent from the equilibrium distributions. 



Zeroth order 



At zeroth order in fluctuations one obtains 
P^a,Q„^(.,rt^-p^^, ^°(-'^^-^--'(-'^\ (15) 



t^fi f \ t^TT Go{x,p) ~ Gcq{x,p) 

p^df,GQ(x,p) -p'^U^ ^ . (16) 



Equations (15) and (161 determine the evolution of 
Qo{x,p) and Go{x,p). 



Instead of solving ilbh and ( 16 ) directly, we take mo- 



ments of these equations. In order to describe (0-f 1)- 
dimensional anisotropic dynamics we take the zeroth and 
first moments of Eqs. (jis]) and ([16]) assuming that the 



(17) 



distributions Q^{x,p) and Go(a;,p]~are given by the co- 
variant version of the Romatschke-Strickland distribution 
[551 [33], namely 

Qo{x,p) = Go{x,p) = fo{x,p), 

where 



fo{x,p) = exp (^-^^{p.ur+ap-vr^ 



Accordingly, we take 



}cq{x,p) = Goq(a;,p) = /eq(a;,p). 



where 



fcq{x,p) = exp 



p-U 

T 



(18) 



(19) 



(20) 



Note that one can also use anisotropic Fermi-Dirac and 
Bose-Einstein distributions for the (anti-)quarks and glu- 
ons, respectively; however, the only change to the final 
result will be the precise value of the isotropic plasma 
frequency of the system. For the sake of simplicity we 
present the case of a Boltzmann distribution and gener- 
alize the results to the full quantum statistical distribu- 
tions in the end. 



The four- vector appearing in ( 18 ) defines the di- 
rection of the beam (z-axis) 



F^=7.(^'., 0,0,1), 7, = (l-«2)-i/2. 



(21) 



We note that the four-vectors C/'' and y satisfy the 
normalization conditions 



C/2 = 1, 



V^ = -l, U-V = 0. 



(22) 



In the local rest frame of the fluid element, and 
have simple forms 



[/^ = (1,0,0,0), = (0,0,0,1). 



(23) 
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For the (0+l)-dimensional boost-invariant expansion 
considered in this paper, we may use 



= (cosh77, 0, 0, sinhr;), 
= (sinh77, O,0,cosh77). 



(24) 



With the assumptions \n\ and (191, the kinetic equa- 



tions (151 and (161 are reduced to a single equation for 



the background distribution 



^^p,u /oi-^P)'U^^P) _ (25) 



1. Zeroth moment of the kinetic equation 



Integrating Eq. ( 25 ) over three-momentum and includ- 

(26) 



ing the internal degrees of freedom we obtain 
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where A'o and iVcq are particle currents Q 



K 



(27) 



A simple calculation performed in the local rest frame 
gives 



no 



.90 



A3 



90_rp3 



(28) 



Here go is the degeneracy factor accounting for internal 
degrees of freedom (we show below that the equations of 
motion for the background are insensitive to the specific 
choice of 5o)- For longitudinal boost-invariant expansion 
one finds 



d 



T 



Thus, using (127|) and (l28j) in (l26j), we obtain 

{T/Af yrre - 1 



3 dA 



1 



1 



AdT 2{1 +C)dT T 



2. First moment of the kinetic equation 



(29) 



(30) 



In the next step we multiply Eq. (25 1 by p'' and inte- 



grate over three-momentum. In this way, we obtain 



(31) 



' eq 



^ There is no term proportional to V in A^^ , due to the quadratic 
dependence of fo on V^. 



where [H [23] 

To"" = (eo + Pi^WU" - P^g^'" - (Fj_ - P\\ W^'V (32) 
and 

= (£cq + PcciWU" ~ P,^g^\ (33) 

In order to conserve energy and momentum, the right- 
hand-side of (31 1 should vanish. Hence, we obtain the 



Landau matching condition 
where 

3.goA4 



350^4 



cq 



and the function 7?.(^) has the form [T5] 



no 



1 



2(1 + 



(1 + OarctanV^ 



(34) 
(35) 

(36) 



Eqs. (34 1 and (35) are used to obtain the ratio T/A 
needed in ( 30 1 



r = A7^l/4(e). 



(37) 



For purely longitudinal boost-invariant motion, the 
energy-momentum conservation law d^^T^'^ — takes a 
simple form 



dep 
dr 



f + P\\ 



(38) 



Eq. ( 38 ) may be reduced to the equation 



7^'(0§ + 47^(e)^ = -I (no + ] , (39) 



where 



no 



(40) 



3. Evolution of the time-evolving background 



Eqs. (30), (37), and (39) provide three equations for 



three unknown functionsri\.(T), ^(r), and T{t). The so- 
lutions of these equations allow us to determine the back- 
ground for the plasma oscillations. Below we will con- 
sider two cases: (i) a fixed relaxation time of r^q = 1 fm/c 
and (ii) a relaxation time that is proportional to the lo- 
cal inverse transverse momentum scale of the plasma. In 
the second case the relaxation time is fixed by requiring 
that, in the limit of small momentum-anisotropy, the lin- 
earized anisotropic hydrodynamics equations reproduce 
2nd-order viscous hydrodynamics [T^. In this case, the 
relaxation time is given by 



' cq 



(r) 



577 



27^l/4(OA ' 



(41) 
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(a) Teq= 1 fm/c 
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(b) AnT] = 1 



drodynamic evolution is taken to be tq = 0.1 fm/c, 
and the final time is taken to be tj — 10 fm/c. 
We consider three different choices of the initial condi- 
tions corresponding to three different values of the ini- 
tial anisotropy — £,[to): = 99 (dashed lines), 
^0 = (sohd lines), and ~ —0.99 (dotted lines). 
The initial transverse-momentum scale Aq is taken to 
be Ao = (1 GeV)-(l -l-^o)^^^ in each case. The factor 
of (1 + ^o)^^^ guarantees that the initial number den- 
sity of the background is the same for all values of 
considered. We show the proper-time evolution of the 
anisotropy parameter, ^, in three different cases: (a) fixed 
Toq = 1 fm/c, (b) Toq givcu by Eq. (41 1 with ^1:7) = 1, and 
(c) Eq. ( [4l| ) with Airf] — 10. In the case of Tcq = 1 fm/c 
shown in Fig. [ija) the anisotropy vanishes at late times. 
When Toq is given by Eq. (41) (see Fig. [T](b,c)) a finite 
anisotropy remains at late times which is consistent with 
2nd-order viscous hydrodynamics. Note that the initial 
growth of the anisotropy is connected with the effects of 
free streaming which dominate the very early dynamics 

mm. 



0.1 - 



B. First order 
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(C) 4^7 = 1 



At first order in fluctuations we obtain 

P^d^5Qf{x,p)±ge, ■ Ff^y^d^Q^ {x,p) 



(42) 



' cq 




p''d^5G,,{x,p) + gri,^ ■ F^>,9PGo(x,p) 
SGij{x,p) 



= -P'^U^ 



(43) 



0.1 - 



In the following we use the boost-invariant variables 
introduced in Refs. |37l |38| 



0.01 W , , 

0.1 0.2 0.5 1.0 2.0 5.0 10.0 

r [fm/cl 



w — tpii — zE, 



and 



(44) 



FIG. 1. (Color online) Time dependence of the anisotropy 
parameter ^ for (a) constant Teq, (b) time varying Tcq with 
i-Kfj = 1, and (c) time varying Toq with Anfi = 10. In each plot 
we show three different initial values of ^: ^0 ~ 99 (dashed 
line), ^0 = (solid line), and ~ —0.99 (dotted line). 



V ^ Et — pn z = \/ w'^ + m'^T^ 



(45) 



From these two equations one can easily find the energy 
and the longitudinal momentum of a particle 



E=p° = 



vt -\- wz 



wt + vz 



(46) 



where fj — rjs/S is the ratio of the shear viscosity to 
entropy density and it is implicitly understood that ^ 
and A depend on proper time. In what follows we will 
assume that fj is time independent. 

In Fig. [1] we show the time dependence of the 
anisotropy parameter ^. The initial time of the hy- 



In addition, we have 

w = Tmj_ sinh {y — rj) , v = Tm± cosh [y — 77) . (47) 

Since the distribution functions are Lorentz scalars, 
they may depend only on r, and w. Therefore, we 



6 



find the general boost-invariant form of the terms ap- 
pearing in Eqs. (42) and (43) 

V 



(48) 



dfo 



The Romatschke-Strickland distribution function takes 
the form 



/o(r, w) = exp 



i-y/pir2 + (l+^(r))«;2 



A(T)r 



Using Eqs. (48) we find the following equations 
— ^Tgrei-S^- 

OT aw Tcq 



(49) 

(50) 
(51) 



The formal solutions of (51 1 for constant or time- 
dependent Tcq may be expressed as integrals 



Here we used the property SQ~ (r, p±,w) = 
6Qf {t,p±,—w) to eliminate the antiquark distribu- 
tion function. 

The invariant measure in momentum space is 

2 '^P\\ 2 

dP ^ Vgd pj_—^ ^ Vgd pj_ — , (58) 

where — t'sf/(27r)'^ and i^sf denotes the number of in- 
ternal degrees of freedom connected with spin or flavor 
{vsf = 4 for quarks and Vgf — 2 for gluons). Using Eqs. ^ 
and ( 58 ) we find 



' J dr'T'Dir^) J d^p^ J ^ 



dwp'^ 



(27r)3?; 



X— (r ,w) 
aw 



2.4. - 5,(r') + 2.3-5,(r') 



For the Lorentz index v = Q, we use the symmetry of 
the distribution function under the change w — > —w and 
write 



5Qf = Tge.- j dr\'D{r,r')£^{r')^{r\w), (52) f = W z 9^ J rfr' Z?(r, r' )£.(/) ^ 



(59) 



6G, 



-gVi 



dT'T'D{T,T')e.iT')^{T',w), (53) 



where we have introduced the damping function 
D[T, T ) = exp 



^cq(r") 



(54) 



V. MAXWELL EQUATIONS 

In order to close our system of equations we have to 
couple the fluctuations via currents to the Maxwell equa- 
tions 

9MF^''=j^ (55) 
where the color current is given by the expression 

» 3 3 

r^g dPp'^ i^Qt - SQi) + E ^^/G^:' 

1=1 i,j = l 

(56) 

Using the formal solutions for the distribution functions 
at first order, we obtain 

T 

f = -g' J dT'T'D{T,T') J dPp'^^iT'.w) (57) 



d^p± f dw w 
(2^ 



+ (1 + 



Here, the primes denote the dependence on r', for exam- 
ple, A' = A(t') and = /o(r',pj^, w). 

The zeroth component of the Maxwell equation gives 



T dr ' 



(60) 



hence, Eqs. ^ and ^ yield 



dr 



(27r)2r 



A' 



To 



OO CO 



X / dpi 



dw w 







« VpI^" + (1 + e) 



w^ 



(61) 



To proceed, we introduce new variables 7 and (f> defined 
by the relations 



7 cos ( 



P± 
A'' 



7sin0=(l + e')^ 



(62) 



3 3 
^2e,(e, . 5,(r')) + J2 ' ^-(^')) 



^ Because of boost-invariance we obtain the same result from the 
Maxwell equations with u = 3. 
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The integration over 7 (from to oo) and (j) (from to 
7r/2) can be performed analytically 1311. and one obtains 



dr T 



T 

J dr' D{t,t')£,{t') 



(63) 



where, for a Boltzmann distribution, k — 10/(37r^) and 
the function J{Ty/l + /t') is defined by the formula [3T] 



36 + 1 

4 63/2 



arcsm 



5-1 



b = 



-12 



(i + O-i- 



3 

26' 
(64) 



Eq. ( 63 ) determines the oscillations of a uniform chro- 



moelectric field in an arbitrary time-evolving anisotropic 
background. It forms the basis of our numerical calcula- 
tions presented below. Interestingly, the form of (631 is 
very similar to that obtained in |31| . Eq. ( 63 ) is reduced 
to Eq. (4.8) in [31^, if we set ^' = in the argument of 
J, A' is replaced by the temperature T', and the damp- 
ing function D is taken to be equal to unity. Note that 
since the calculations in j31j were performed using the 
quantum statistical distribution functions, the factor k, 
in Ref. [H] equals [N^ + iV//2)/9 = 4/9 for JV^ = 3 and 
Nf = 2. The calculations presented in this paper can 
be also performed using quantum statistical distribution 
functions and, in that case, one obtains k, = 4/9. The 
results presented in Sec. IVlluse this choice of n. 



A. Special Case: Longitudinal free streaming limit 

For the case of longitudinal free streaming it is possi- 
ble to write the integro-differential equation (63) as an 

00 



ordinary differential equation. In the limit T( 



one has ^(r) 



{T/n^oY " 1 and A(t) 



cq 



Ao with 

^(tq) being the point in time when ^ = 0. 



Tiso = To/ V 1 

Inserting these relations into ( 63 1 gives 



J 



dr 



r 

J dr'e.ir'y, 



(65) 



where = ng^A^. Taking a derivative of both sides of 
this equation with respect to r gives an ordinary differ- 
ential equation 




dS^ir) 



dr 



= -^1,S.{t). (66) 



which is supplemented by the initial conditions £z{to) — 
Szfl and £'^{tq) = 0, where the latter condition follows 
from (65) upon setting r = tq. 



The differential equation above is nonlinear and must 
be solved numerically; however, at late times we can find 
the asymptotic form of the solution by using 



lim J 

>oo V Tip 



to obtain 



1 d 
T dr 



At 



idSzir) 
dr 



O 



-M^2(t) 



(67) 



(68) 



where /i — 37rTiso'^pi/4. This differential equation has a 
solution of the form 



lim £, 



1 



AJi {2y/JI¥) + BYi {2y/JI¥)] , (69) 



where Ji and Yi are Bessel functions of the first and 
second kind, respectively, and A and B collect undeter- 
mined constants that will be matched below. Note that 
this agrees with the result first obtained in App. A sub- 
section 2 of Ref. [31] 



VI. RESULTS 

For the numerical results we choose a particular direc- 
tion of the chromolectric field by aligning it initially in 
the '3' direction, i.e. £z{tq) — {Sz (to),0). The evolu- 
tion keeps £z in the '3' direction at all times. In Fig. [2] 
we plot the time dependence of the normalized longitu- 
dinal chromoelectric field £i for (a) constant r^q, (b) 
time varying Toq with Anfj = 1, and (c) time varying 
Toq with Anfj = 10. We turn off the damping manually 
by setting the damping function D{t,t') = 1. In each 
plot we show three different initial values of ^: ^0 = 99 
(dashed line), Co = (sohd hne), and = —0.99 (dotted 
line). For each we have adjusted Ag as described in 
Sec. |IV A3 in order to guarantee that the initial num- 
ber density is held constant. This figure demonstrates 
that, although we have varied our assumed value of ^g 
over a large range corresponding to initially extremely 
prolate (^g = —0.99) to extremely oblate (^g = 99), if 
the results are normalized such that the initial number 
densities are held constant, the resulting oscillations are 
not dramatically different. This is due to the fact that 
what sets the time scale for the plasma oscillations is 
Wpi and, generally, one has Wpi ^ n(rg)/Ag. However, 
despite being qualitatively similar, there are important 
quantitative differences which remain. 

In Fig. |3] we plot the time dependence of the normal- 
ized longitudinal chromoelectric field £z for (a) constant 



^ Their result is expressed in terms of the longitudinal vector po- 
tential. One must compute the longitudinal electric field using 
= n'' = T^^drAn in order to compare with our result. 
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FIG. 2. (Color online) Time dependence normalized of the 
normalized longitudinal chromoelectric field £z for (a) con- 
stant Teq, (b) time varying Tcq with i-Kfj — 1, and (c) time 
varying Tcq with 4nf] = 10. We turn of! the damping manually 
by setting the damping function D{t,t') = 1. In each plot 
we show three different initial values of ^: ~ 99 (dashed 
line), ^0 = (solid line), and ~ —0.99 (dotted line). 



FIG. 3. (Color online) Time dependence of the normalized 
longitudinal chromoelectric field £z for (a) constant Teq, (b) 
time varying r^q with infj = 1, and (c) time varying Teq with 
4:TTfj = 10. In each plot we show three different initial values 
of C: Co = 99 (dashed line), = (solid fine), and Co = -0.99 
(dotted line). 



Toq, (b) time varying Teq with Anfj = 1, and (c) time 
varying Toq with ^-kt] = 10. We now include the damp- 
ing function in the integrand of the integro-difFerential 
equation. In each plot we show three different initial 
values of ^: = 99 (dashed line), = (solid line). 



and ^0 = —0.99 (dotted Une). Once again, for each 
we have adjusted Ao as described in Sec. |IV A3| in or- 
der to guarantee that the initial number density is held 
constant. From this figure we see that the damping func- 
tion D{t,t') has an extremely important impact on the 
time evolution of the oscillations of uniform longitudinal 
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FIG. 4. (Color online) Comparison of solution to the ordi- 
nary differential equation ( 66 1 (thick black line) , the integro- 



differential equation (63 1 (thin red line), and the asymptotic 



solution ( 69 1 (thick gray dashed line) in the case of longitu- 
dinal free streaming (rcq = oo). The inset shows the three 
results scaled by r'^^^ . 



chromoelectric fields|f] 

In Fig. |4] we compare solutions to the ordinary dif- 
ferential equation ( |66|) (thick black line), the integro- 
differential equation ( |63[ ) (thin red line), and the asymp- 
totic solution ( 69 ) (thick gray dashed line) in the case of 
longitudinal free streaming (req = oo). The inset shows 
the three results scaled by r'l'^ in order to make a more 
precise comparison of the late time behavior. For the 
asymptotic solution the unknown constants appearing in 
( 69 1 were fixed by matching numerically to the full solu- 
tion at T = 10'* fm/c. 



As can be seen from this figure, 
our numerical solution to the integro-differential equa- 
tion ( 63 1 gives the same result as direct solution of the 



differential equation ( 66 ) and at late times both agree 
well with the asymptotic solution. This gives us confi- 
dence that the numerical method used for solution of the 



integro-differential equation (63 1, in the general case, is 
reliable. 



VII. CONCLUSIONS 

In this paper we studied the oscillations of a uni- 
form longitudinal chromoelectric field in a dynamically- 
evolving momentum-space anisotropic background. The 
required anisotropic hydrodynamics equations were ob- 
tained from the first two moments on the Boltzmann- 
Vlasov equation using a spheroidal form for the local 



rest frame one-particle distribution function. The result- 
ing anisotropic hydrodynamics equations provided the 
proper-time dependence of the local transverse momen- 
tum scale, A(t), and momentum-space anisotropy, ^(t). 
We then expanded the Boltzmann-Vlasov equation to 
first order in fluctuations and coupled these fluctuations 
to the Maxwell equations. From this procedure we ob- 
tained an integro-differential equation which governs the 
time evolution of a uniform longitudinal chromoelectric 
fleld. The integro-differential equation allows for an arbi- 
trary time dependence of the scale A(r) and momentum- 
space anisotropy ^(t) provided that the system is boost 
invariant and transversely homogeneous. The integro- 
differential equation obtained also includes the effect of 
coUisional damping of the oscillation. 

Having obtained the integro-differential equation nec- 
essary, we proceeded to solve it numerically for a vari- 
ety of different initial momentum-space anisotropics us- 
ing two different assumptions for the relaxation time Toq. 
We considered (i) the case of constant Toq, in which case 
the late-time dynamics is that of ideal hydrodynamics, 
and (ii) a time-dependent Teq that is inversely propor- 
tional to the local inverse average transverse momen- 
tum, in which case the late-time dynamics is that of 
2nd-order viscous hydrodynamics. We showed that for 
flxed initial number density the effect of time-varying 
momentum-space anisotropy is important but not over- 
whelming large. However, we found the effect of coUi- 
sional damping to be quite important for the dynamics 
of stable chromoelectric oscillations. 

We should stress that the results presented here are 
exploratory in the sense that we have only studied the 
dynamics of a stable uniform chromoelectric field; how- 
ever, the general method used here could have a wide- 
ranging application in the study of the dynamics of all 
stable and unstable modes of a longitudinally expand- 
ing QGP. Previous studies in this direction have been 
restricted to the limiting case of a longitudinally free- 
streaming anisotropic background [32H31]- We were able 
to check our results against the free-streaming results ob- 
tained by Rebhan and Steineder [34] and found that we 
are able to reproduce their results for the dynamics of 
a uniform longitudinal chromoelectric field. In addition, 
we were able to express the free-streaming evolution of 
a uniform chromoelectric field as an ordinary differential 
equation, albeit a highly-nonlinear one. It would be very 
interesting to see if the results obtained here could be 
used to obtain modified hard-loop equations of motion 
that would allow one to simulate the full non-Abelian 
dynamics of stable and unstable modes in a realistically 
time-evolving anisotropic background. Work along these 
lines is in progress. 



* We made a preliminary study of the impact of coUisional damp- 
ing on unstable modes and found that the damping serves only 
to slightly reduce the growth rate of unstable modes. Results 
from this study will be reported elsewhere. 
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